Searching for Spurious Solar and Sky Lines in the Fermi-LAT Spectrum 
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We search for a unified instrumental explanation of the spectral features seen near E 1 — 130 
GeV in photons collected by Fermi-LAT from the galactic center and from the Earth's limb. We 
report for the first time a similar feature in photons originating from the vicinity of the Sun, and 
examine the instrumental characteristics of this Solar feature. To test an instrumental hypothesis, 
we identify the range of photon incident angles where most of the peak photons are observed in 
these three spectral features. An examination of the spectrum of photons from the rest of the sky 
with this characteristic angular range reveals a hint of a spectral feature near E-, — 130 GeV. These 
results cast further doubt on the dark-matter-annihilation interpretation of the galactic center peak. 



INTRODUCTION 

The particle nature of dark matter remains a mystery. 
One potential avenue for discovery is via dark matter 
annihilation into standard model particles. If the an- 
nihilation results in fermions or heavy bosons, photons 
may be produced via hadronization and the decay of 7r° 
particles. These photons are expected to be fairly low- 
energy {E~f <ps 50 GeV) and are expected to be diffi- 
cult to distinguish from other sources. A clearer feature 
may appear from annhilation directly into two-body final 
states including a photon. Rather than yielding a broad 
energy spectrum, this process would produce a photon 
with a well-defined energy simply related to the mass 
of the dark matter particle. Such high-energy gamma 
rays typically do not scatter in transit to the Earth from 
the regions of high dark-matter density where they are 
produced, making their energy and direction useful han- 
dles. The Fermi Large Area Telescope (LAT) has a large 
field-of-view and excellent energy resolution over a broad 
range of energy (20 MeV - 1 TeV) Q, making it a pow- 
erful probe for peaks in the photon spectrum from dark 
matter annihilation 0, 

Recently, a statistically significant peak has been re- 
ported in the Fermi-LAT photon spectrum near E 1 = 130 
GeV [3, Q with a source close to the galactic center 0- 
Q. The importance of such a discovery, if interpretted 
as due to dark matter annihilation, requires a careful 
exploration of other more mundane explanations, such 
as potential astrophysical spectral features in the non- 
dark-matter background from the galactic center, or in- 
strumental effects in the Fermi-LAT detector. Suspicious 
features have been reported in the incident angle of the 
photons from the galactic center Q and a statistically 
significant feature has been reported in photons from the 
Earth's limb [loj ]. where no signal is expected from dark 
matter annihilation. 

In this paper, we consider another luminous source, 
the Sun, for evidence of spectral features. In addition, 
we use the sPlots lli] algorithm to reconstruct - sepa- 
rately for peak and background photons - the distribu- 



tions of instrumental quantities in an attempt to iden- 
tify a common instrumental characteristic unique to the 
observed spectral features. Such a common characteris- 
tic would point strongly to an instrumental explanation 
which would be independent of the source of the pho- 
tons. To test such potential instrumental explanations, 
we examine the spectrum of photons with these instru- 
mental characteristics from the entire sky for evidence of 
spurious features. 

THE FERMI-LAT DATA SAMPLE 

We use the publically available Fermi-LAT photon 
data through Jan uary 3rd 2013, making standard quality 
requirements fl2l . fl3l | . In addition, we consider several 
additional requirements for the distinct samples. The 
analysis of the Earth's limb photons requires a zenith 
angle (measured with respect to the zenith line, which 
passed through the earth and the Fermi spacecraft) be- 
tween 110° and 114° degrees; the analysis of the Solar 
photons requires a zenith angle less than 105°, rocking 
angle of less than 52°, and angular distance from the po- 
sition of the Sun of less than five degrees. The galactic 
center spectrum is drawn from a five-degree circle around 
(I = 0, b = 0) in galactic coordinates, requiring zenith an- 
gle less than 105° and rocking angle of less than 52°. 

Other than the reconstructed energy, the photons have 
other measured characteristics [3| which may give in- 
sight into instrumental effects: 

• incident angle 9, measured with respect to the top- 
face normal of the LAT, 

• azimuth angle 0, measured with respect to the top- 
face normal of the LAT, folded as described in Eq. 
(15) of Ref. 0. 

• mission elapsed time, measured relative to January 
1, 2001, 

• conversion type (front or back), indicates whether 
the event induced pair production in the front 
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(thin) layers or the back (thick) layers of the 
tracker, 

• the magnetic field in which the LAT is immersed, 
as parameterized by the Mcllwain B and L param- 
eters |15| . 



SPECTRA 

For the observed feature we assume a single line where 
the pdf /n n e (E 7 |i?iine) is defined according to the Fermi- 
LAT energy dispersion tools definition 116] with a true 
photon energy of En ne (see Fig. 3 of Rcf [9j). To analyze 
the energy dependence, we integrate out the dependence 
on 9 using the exposure and livetime in the region of 
interest. 

The background pdf is a simple power-law: 



/ bg (£ 7 |/?,a) =/3 
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FIG. 1: Observed Fermi-LAT photon spectrum from the 
Earth's limb, requiring in addition 30° < p hoton < 45°, fol- 
lowing Ref. [J. 



Earth's Limb Spectrum 

Some fraction of the time, the LAT's field of view in- 
cludes the limb of the Earth, where photons are produced 
in collisions of cosmic rays with the atmosphere. These 
are not salient for dark matter searches, given their ter- 
restrial source, and so are typically discarded. The obser- 
vation of a statistically significant feature near E~ = 130 
GeV in this spectra for photons with 30° < 6 < 45° 
(Ref [7( and Fig. [1]) where none is expected from dark 
matter annhiliation points strongly to an instrumental 
explanation of this feature. 

Solar Spectrum 

In the same spirit, we examine the high-energy So- 
lar spectrum. The Fermi collaboration studied the low- 
energy (J5 7 < 10 GeV) spectrum in a smaller dataset fl7j ]. 

While arguments have been made that the Sun's grav- 
itational well may trap dark matter particles [HI], and 
searches have been performed for evidence of Solar dark 
matter annhiliation via neutrinos (l9j . the majority of 
the dark matter is expected to be gathered near the cen- 
ter of the Sun. Photons resulting from Solar dark-matter 
annhiliation would not escape the Sun without interact- 
ing, and so no signal is expected in the Solar photon 
spectrum [20[ unless the Solar dark matter halo is much 
broader, contrary to current ideas. A feature in the Solar 
spectrum would suggest instrumental causes, in the same 
manner as the peak in the Earth's limb data. 

To isolate Solar photons, we use the position of the Sun 
relative to the LAT and calculate the effective position 



of the Sun in galactic coordinates at the moment a given 
photon is received, (&q,Zq). The Solar region is then 
defined via the angular distance 

cos(A_Rq) = sin(6 7 ) sin(&Q)+cos(& 7 ) cos(&q) cos(Z 7 — Zq) 

We consider photons with AR@ < 5° outside the re- 
gion of the galactic center (see Fig. [3J left) to be Solar 
photons. The Solar photon spectrum is shown in Fig. [2] 
(right), compared to the spectrum of the inverted region 
(ARq > 5°, same veto on galactic center photons). The 
local statistical significance of the feature near E 1 = 135 
GeV is 3.2er (see Fig. 2]) which does not suffer from a 
look-elsewhere effect, as the region-of-interest in E 1 of 
the feature is determined by the feature in the galactic 
center spectrum. This appears to be the first examina- 
tion of the high-energy Solar spectrum {E 1 > 10 GeV) 
using Fermi-LAT data. The same feature appears for 
different requirements of ARq, see Fig. [3] in fact, mask- 
ing out the physical extent of the Sun (ARq € [1, 5]°) 
yields nearly the same spectrum, suggesting that these 
photons do not originate from the Sun proper, but from 
its vicinity. 

INSTRUMENTAL STUDY WITH SPLOTS 

To consider potential instrumental explanations of the 
observed features, we ask whether the photons in the 
peaks have distinct instrumental characteristics from the 
background photons. A naive approach would be to sim- 
ply divide the photons into two categories by energy 
range; this leaves the peak-region sample with signifi- 
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FIG. 2: Left, galactic coordinates of photons within ARq < 
5° of the position of the Sun. Right, energy spectrum for 
photons from the vicinity of the Sun (A_Rq < 5°, excluding 
the galactic center), compared to photons from the rest of the 
sky (ARq > 5°, excluding the galactic center), normalized 
to equal yield. 
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FIG. 3: Clockwise from top left: distribution of ARq for 
photons with _B 7 6 [125, 140] GeV; energy spectra for photons 
from the vicinity of the Sun using ARq < 3° , ARq < 10° , 
or ARq £ [1, 5]°, excluding the galactic center, compared to 
photons from the rest of the sky, again excluding the galactic 
center. 



FIG. 4: Energy spectrum as shown in Fig.[5]for photons from 
the vicinity of the Sun (ARq < 5°), excluding the galac- 
tic center, with background-only and background-plus-signal 
hypotheses. 
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FIG. 5: Energy of Fermi-LAT photons with signal plus back- 
ground fit. Left is the spectrum from the Earth's limb with 
9 < 60° , right is the Solar spectrum with ARq < 5° , exclud- 
ing the galactic center. 

data and application to the Fermi-LAT galactic center 
spectrum, see Ref [9j. 



cant background contamination. A superior technique is 
statistical unfolding using the sPlots algorithm [llj . 

If a sample of events contains contributions from mul- 
tiple sources, whose relative contributions can be mea- 
sured by fitting pdfs in a discriminating variable, the 
sPlots algorithm can statistically extract the distribu- 
tion of each of the sources in other variables of interest, 
which we refer to as the 'unfolding variables'. Note that 
sPlots requires knowledge of the pdf in the discriminat- 
ing variable, but no assumptions are required of the pdfs 
in the unfolding variable other than that the pdfs can be 
factorized between the discriminating and unfolding vari- 
ables. For a brief derivation of SPlots, examples in toy 



Analysis 

To analyze the features of the Fermi-LAT data using 
sPlots, we must fit the background and signal pdfs de- 
scribed above in the discriminating variable, £L. Ap- 
plying these pdfs to the observed limb and Solar photon 
energy spectrum yields the fits seen in Fig. [5] In the 
limb spectrum, we have not required 30° < 9 < 45° as is 
done in Fig.[TJ but instead we apply a looser requirement 
6 < 60° in order to study the (9-dependence. 

Unfolded distributions of incidence angles are shown 
in Fig. [6l The recorded time and conversion type are 
in Fig. [7J The magnetic field parameters are shown in 
Fig. 1 
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FIG. 6: Disentangled signal and background distributions. 
Top, cos(f3) where 9 is the photon incidence angle relative to 
a line normal the Fermi-LAT face. Bottom, <j>, the photon 
incidence angle relative to the sun- facing side [1J]. Left are 
photons from the Earth's limb, right are photons from the 
vicinity of the Sun. For the Solar spectrum, 4> is aligned with 
respect to Fermi's solar panels, which track the Sun. 

In each case, we compare the distributions quantita- 
tively by calculating the x 2 /dof between the peak and 
background distributions. As the signal and background 
weights are anti-correlated, this is calculated as 

2= V (^cak"^g.) 2 

where N^ cak is the sum of the weights ^ sP poa k in that 
bin, and A7Vp eak is calculated from toy simulations which 
estimate the expected variance of the measurement in 
each bin. Similar expressions apply for the background 
uncertainties. 



Discussion 

Our previous analysis of the galactic center feature 
showed that a large fraction of peak photons had cos(6*) 
near 0.7, unlike the photons from the background (Fig. 
5 in Ref. [9(). The limb feature photons also peak near 
this region, first reported in Ref [7], and confirmed here, 
see Fig. |5l The Solar spectrum statistics are too poor to 
contribute to this question. 

The galactic center instrumental analysis also showed 
some discrepancy in the magnetic field environement of 
the LAT when the photons were recorded (Fig. 11 in 



FIG. 7: Disentangled signal and background distributions. 
Top, the mission elapsed time since Jan 1 2001 [H|] . Bottom, 
fraction of events in which the pair production is induced in 
the front (thin) or back (thick) layers of the tracker. Left are 
photons from the Earth's limb, right are photons from the 
vicinity of the Sun. 

Ref. @), near 1.6 Gauss. The limb photons have a mi- 
nor discrepancy near 1.5 Gauss, but entirely consistent 
with statistical fluctuations. Similarly, the galactic cen- 
ter photons have a feature near folded azimuth <f> = 0.8; 
the limb photons have some features near <fi — 0.7 and 
<j> = 0.9, but a deficit near </> = 0.8 with respect to the 
background photons. 

THE REST OF THE SKY 

If the feature near E 1 = 130 GeV has an instrumental 
rather than astrophysical explanation, then it should be 
independent of the source of the photons and one should 
be able to identify similar features in any spectrum. The 
full-sky spectrum with no restrictions on instrumental 
characteristics does not show such a feature; if it exists, 
a feature in the full-sky spectrum may be localized to 
particular instrumental regions and washed out in the 
unrestricted spectrum. 

To discover potential instrumental causes, we look for 
universal characteristics of the peaks found in the three 
spectra. If the full sky spectrum in a restricted range of 
these instrumental characteristics were to reveal a fea- 
ture near E 1 = 130 GeV, it would lend support to the 
instrumental explanation. 

We first consider the incident angle, 8, already shown 
to reveal the Earth's limb peak. In Fig. [9j we show the 
spectra in the three regions defined earlier: the galac- 
tic center, the Sun, and the Earth's limb. In addition, 
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FIG. 8: Disentangled signal and background distributions. 
Top, magnetic field strength in terms of the Mcllwain B pa- 
rameter. Bottom, the Mcllwain L parameter. Left are pho- 
tons from the Earth's limb, right are photons from the vicinity 
of the Sun. 



we show the spectrum in the rest of the sky, defined 
by requiring zenith angle < 105°, rocking angle < 52°, 
Ai?Q > 5°, and \Jb 2 + 1 2 > 5°. These four regions are 
disjoint. 

First, we examine the spectrum with no restrictions 
on 9. There are features in the galactic center and Solar 
spectrum, but nothing evident in the limb or sky. When 
we restrict 9 to [30,45]°, the limb feature is evident, and 
we appear to capture a large fraction of the galactic cen- 
ter and Solar features. Most interestingly, a small feature 
is evident in the sky spectrum near E 1 = 130 GeV. This 
feature is intriguing, but is not large enough to unam- 
biguously confirm an instrumental issue in the region 9 
to [30,45]°. 

A similar analysis using the detector azimuth (Fig. [TO]) 
is able to capture most of the galactic center and Solar 
spectra, but do not reveal features in the limb or sky 
spectra. 



not dark matter annihiliation. 

We also analyze the instrumental characteristics of the 
Earth's limb feature, and attempt to identify a universal 
characteristic of photons in the three observed features 
using the sPlots algorithm. 

Finally, we show that a narrow range of incident pho- 
ton angles includes many photons from the observed 
peaks, and reveals a feature near E 7 = 130 GeV in the 
spectrum of the remainder of the sky. While we attempt 
to remove dependence on photon source by examining 
the full sky and placing restrictions only on 9, instru- 
mental characteristics may in general be correlated with 
the photon source, as the geometry of the Sun, Earth and 
galactic center relative to the LAT may give different av- 
erage distributions of instrumental characteristics. If the 
spectral feature were caused by an instrumental issue in 
a slice of 9 and another not-yet-identified variable, it is 
plausible that photons from the Earth, Sun and galactic 
center sweep across this slice differently, giving different 
^-dependences and explaining the stronger ^-dependence 
of the limb and sky features than in the solar and galactic 
center features. In addition, more precisely focusing on 
photons in such a slice would enhance the feature in the 
sky spectrum. 

While this does not give a definitive instrumental ex- 
planation for these spectral features, it casts significant 
further doubts on the dark matter hypothesis. 
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CONCLUSIONS 

We report a spectral feature near E 1 = 130 GeV in 
the spectrum of photons from the vicinity of the Sun, 
where none is expected from dark matter annhiliations. 
The statistical significance is 3.2c, corresponding to a 
p-value of < 0.5% that the feature is a statistical fluctu- 
ation. This is the first report of such a feature. Based on 
current understanding of the Solar dark matter halo, if 
the feature in the Solar spectrum is more than statistical 
fluctuation, it represents evidence of instrumental issues, 
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FIG. 9: Energy spectra for four regions; from left: galactic center, Solar vicinity, Earth's limb, and the full sky with the GC 
and Sun removed. Top shows the complete spectrum, center is for photons with incident angle 6 in the range [30, 45] degrees, 
and bottom for photons with 6 out of the range [30, 45] degrees. 
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FIG. 10: Energy spectra for four regions; from left: galactic center, Solar vicinity, Earth's limb, and the full sky with the GC 
and Sun removed. Top is for photons with azimuth angle (f) in the range [—160, 50] degrees, and bottom for photons with (j> out 
of the range [—160,50] degrees. 
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